Linear and nonlinear wave propagation in weakly relativistic quantum plasmas 
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We consider a recently derived kinetic model for weakly relativistic quantum plasmas. We find 
that that the effects of spin-orbit interaction and Thomas precession may alter the linear disper- 
sion relation for a magnetized plasma in case of high plasma densities and/or strong magnetic fields. 
Furthermore, the ponderomotive force induced by an electromagnetic pulse is studied for an unmag- 
netized plasma. It turns out that for this case the spin-orbit interaction always give a significant 
contribution to the quantum part of the ponderomotive force. 
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I. INTRODUCTION 



Much recent work have studied quantum effects in plasmas. Applications typically occur for high density plasmas 
of low or moderate temperature, see e.g. Refs [3-[l| for detailed discussions. Typically hydrodynamic approaches 
cover effects due to particle dispersion and the Fermi pressure @, § > whereas much of the kinetic treatments is based 
on the Wigner equation [l[ . More accurate treatments based on the Kadanoff-Baym kinetic equations p| Q are also 
common. Moreover, the magnetization currents and magnetic dipole force due to the electron spin have been included 
in hydrodynamic Q as well as kinetic theories [llj- 

Relativistic effects may also have considerable influence over the dynamics due to e.g high power lasers (l2^ . 
Combined quantum mechanical and relativistic effects in plasmas might be possible to probe with soon to be built 
lasers, in particular x-ray lasers like X-FEL (e.g [l3j]) and future lasers with even higher photon energies. Also in 
astrophysics there are domains where quantum relativistic plasma dynamics are of importance [Tij . 

In the present paper we will study wave propagation using a recently developed kinetic model [l5j . extending 
previous quantum mechanical models (To| . to including weakly relativistic effects such as spin-orbit interaction and 
Thomas precession [l6j]. Firstly the dispersion relation for electromagnetic waves propagating parallell to an external 
magnetic field is derived. It is found that the new terms modify the dispersion relation for high densities and/or 
strong magnetic fields. The linear investigations is then taken as a starting point for calculating the ponderomotive 
force due to electromagnetic pulses. The classical version of this problem has been thoroughly studied in the past 
decades [l7|, and the ponderomitive force is known to give raise to effects such as wakefield generation [l8[, soliton 
formation, self-focusing and wave collapse 1 1 ll| . The contributions to the ponderomotive force from non-relativistic 
spin dynamics have been studied by Refs. [2l|, [22|], and the it has been shown that the this may lead to spin- 
polarization of an initially unpolarized plasma. For the specific case of an unmagnetized plasma we find that the 
weakly relativistic effects significantly modify the spin contribution to the ponderomotive force. Furthemore, the spin 
quantum contribution becomes comparable to the classical ponderomitive force when the electromagnetic wavelength 
approaches the Compton wavelength. The implications of are results are discussed. 



II. BASIC MODEL AND LINEAR SOLUTION 



In a recent work Asenjo et al |15| presented a kinetic evolution equation for a weakly relativistic spin 1/2 collisionless 
plasma in the long scale lenght limit, to order c~ 2 and h 2 : 
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where / = j (x,p, s,t) is the quasi-distribution function defined on a phase space extended by two spin dimensions 
(denoted s) on the unit sphere, in addition to the traditional space and momentum coordinates x and p. We use the 
notation m for the mass and p, is the magnetic moment of the particle, q is its charge. The index x, p or s on the 
nabla operator indicates that it acts on the respective coordinates. When writing down this equation the last term of 
the equation derived in [l5[ was omitted. This term is associated with the Darwin term in the Hamiltonian, and was 
dropped since in a long wavelength expansion it is smaller than the other terms. 
This Vlasov-like equation is coupled to Maxwells equations 
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where the total charge and current density are given by 
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Here pp = qjdflf is the free charge density and the free current density, the polarisation and magnetisation are 
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We make the division / = /o + / where /o is the background distribution and / is the perturbed distribution function, 
assumed to be homogenious and isotropic in the momentum variable. To sum up our model, it is a Vlasov-like equation 
for a quasi distribution function in a phase space expanded by the spin variable s, which is measured in the rest frame 
of the particle. It contains the Lorentz force, magnetic dipole force, Thomas correction and spin precession with spin 
orbit correction. It should be noted that since this model is semirelativistic, the relation between the momentum p 
and the kinetic velocity v is nontrivial. 
After linearization, ([I]) reads 
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We use a standard ansatz of quasi- monochromatic harmonic variation on the perturbed quantities, E\ = E\ expi(k ■ 
x — ujt) etc, and choose k = kz and the polarization in the (x, y)-plane. The external magnetic field Bq is assumed 
to be static, homogeneous and point in the z-direction. For the momentum variable we use cylindrical coordinates, 
and for the spin we use spherical coordinates on the unit sphere. Furthermore the unperturbed spins are assumed to 
be in thermal equilibrium, and thus fo(s,p) = fo(p) [1 + tanh (pBo/ksT) cos9 s ] f23| . 
Now we expand f\ in eigenfunctions 
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where c.c. stands for complex conjugate and insert (|10[) in Q. Mutiplying with 1 / (2n) exp[i(n(f> p + m(f> s )] and 
integrating both these angles from to 2tt we can solve for W a ,b- Plugging this back into (fTOf we get the disturbed 
distribution function to first order in perturbed quantities as 
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where uj c = qBo/mc is the cyclotron frequency, uj cg = (g/2)u> c is the spin precession frequency and E± = E x ± iE y . 
In combination with Maxwell's equations and the expressions for the currents, we obtain the dispersion relation 
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As expected the dispersion relation has two solutions, corresponding to left and right circular polarisation. 
Taking the long wavelenght limit k — > the result simplifies to 
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Here we can note that for high densities and low temperatures the terms stemming from relativistic effects can 
significantly alter the dispersion relation compared to previous nonrelativistic results like e.g. with [l9| . However we 
will now mainly be concerned with the case where the impact of the new terms are small in a linear calculation, but 
might be of importance when dealing with nonlinear problems. 



III. THE PONDEROMOTIVE FORCE 



The classical pondermotive force has been thoroughly studied, and recently pure spin effects have also been explored. 
In the present paper we are mostly concerned with the effects arising from spin-orbit coupling, which can be considered 
as a quantum relativistic effect. To see these effects, it suffices to study an unmagnetised plasma, thus also reducing 
the algebra to more manageable proportions. This also allows us to make the very reasonable assumption that the 
equilibrium distribution function does not depend on 6 S . This can be motivated since we do not have any external 
field, so the spin states are degenerate. This reduces the first order distribution function to 
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Now we study the evolution equation to second order in perturbed quantities, and only keep source terms on the low 
frequency time scale thus obtaining 
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where the star denotes complex conjugate and the index If indicates that the quantity has no rapid oscillations, and 
thus varies only on the slow spatial and temporal scales. What we want to calculate is the total low frequency current, 
composed in our case by the free current 
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and the polarisation current 
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We note that the magnetisation current vanishes in our geometry. For simplicity we only consider the quantum 
contributions, since the classical ponderomotive current has been calculated in a number of previous works already. 
Starting with the free current we note that 
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Using (|16j) in combination with (|15jl the we can calculate the first term on the right hand side in terms of the field. 
The integral on the left hand side can be dealt with by noting that the convective derivative in the evolution equation 
for the low frequency distribution function is small in the low temperature limit, and can thus be calculated using 
perturbation theory. Now we obtain an expression for the second time derivative of the free current: 
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The approximation performed to obtain (|20[) is the addition of the last term proportional to (p z /m)d/dz. This 
addition is a higher order thermal correction, but is useful since it enable us to rewrite the terms involving fif by 
using Eq (|16[) combined with (|15|) to obtain a driving term for the low-frequency current proportional to the high- 
frequency wave intensity (i.e. proportional to (£^±1 = \E X \ + \E y \ ) Following the same approximate procedure for 
the polarisation current we have 
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Using again eqs (JTHJ) and (fT5|) to to rewrite the right hand source terms and combining the results from (|2U|) and (|2T|) 
we find that the second order time derivative for the total current is given by 
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where we have expanded each term to lowest order in p z to be able to perform the integration (which is consistent with 
the approximations in (|20p and ([21]) ), and defined the plasma frequency lu 2 = A-Kq 2 n^/m. The classical contribution 
has been omitted for simplicity. Using the time derivative of Ampere's law we obtain 
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We can note that the second to last term in the square bracket is what was obtained in previous works based on 
models not containing relativistic effects j22|. If we assume that kc/uj is roughly of order 1 we see that all terms in 
the square brackets are of the same order. This implies that when dealing with an unmagnetised plasma where spin 
effects are important, the spin orbit coupling contributions must be taken into account as well. Furthermore, in the 
approximation where uj 3> uj p such that d/dt = cd/dz we have 
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This spin contribution should be compared with the classical current given by 

and we see that these two source terms will be comparable if hk/mc ~ 1. For typical parameters where u> ~ kc this 
implies that we need photon energies of the order of the electron rest mass energy, i.e. gamma rays. Here it should 
be stressed that for such short wavelengths several other effects that have been omitted is likely to be of importance, 
for example particle dispersive effects and the Darwin term [l5| . 



IV. SUMMARY AND CONCLUSIONS 



In this work we have first solved the linear problem and found the dispersion relation for waves propagating parallel 
to the external magnetic field in a weakly relativistic spin plasma, and seen that it gives the correct classical limit. 
It was also seen that in the long wavelength limit the dispersion relation is significantly modified by the spin-orbit 
interaction when the Zeeman energy approaches the electron rest mass energy. 

Furthermore the nonlinear ponderomotive force in an unmagnetized setting has been derived, and compared with 
previous classical and nonrelativistic quantum result. In the case of high energy radiation the quantum contributions 
are seen to actually dominate over the classical ones. For example if a plasma is illuminated by gamma rays in the 
MeV-regime. The planned high power x-ray lasers like Xfel promise photon energies around 25keV, [l3| which is still 
not enough for the quantum terms to dominate. However, accelerating an electron bunch with a lineac might suffice 
to blue-shift the x-ray photons enough to give them MeV-energies in the rest frame of the electrons. 

Another possibility to obtain the energetic photons is from Gamma ray bursts , see eg (24j for a review. In the 
GRB itself the radiation is created in a relativistic jet, which means that the photons are redshifted in the reference 
frame of the plasma. However there is a possibility that the radiation passes through an accretion disc between the 
source and the observer, and when this happens the ponderomotive forces and subsequent acceleration of the particles 
in the accretion disc can be dominated by the terms calculated above. 

When dealing with highly energetic photons we always face the possibility of QED effects, which are not included 
in the current model. To start with we consider Compton scattering. Studying the Klein-Nishina cross section (25j 
we see that the cross section actually decreases with the photon energy for energetic photons, thus this mechanism 
will be suppressed. Pair production on the other hand can play a dominating part if photon energies are high, but 
in order to conserve energy and momentum the two interacting photons need to be of different energy. If we only 
consider quasi-monochromatic beams or beams with a narrow energy spectrum the photons will therefore not produce 
pairs, and we conclude that it is consistent in this case to neglect QED effects. Furthermore, the particle dispersive 
effects neglected in the kinetic model may also play a role for high photon energies. 
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